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Using discrete element methods, we study numerically the 
dynamics of the size segregation process of binary particle 
mixtures in three-dimensional rotating drums, operated in the 
continuous flow regime. Particle rotations are included and we 
focus on different volume filling fractions of the drum to study 
the interplay between the competing phenomena of mixing 
and segregation. It is found that segregation is best for a 
more than half-filled drum due to the non-zero width of the 
fluidized layer. For different particle size ratios, it is found 
that radial segregation occurs for any arbitrary small particle 
size difference and the final amount of segregation shows a 
linear dependence on the size ratio of the two particle species. 
To quantify the interplay between segregation and mixing, we 
investigate the dynamics of the center of mass positions for 
each particle component. Starting with initially separated 
particle groups we find that no mixing of the component is 
necessary in order to obtain a radially segregated core. 
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I. INTRODUCTION 

When granular materials are placed in rotating cylin- 
ders, different flow dynamics are observed. The ma- 
jor portion of the particles undergoes a solid body ro- 
tation by following the cylinder motion. Close to the 
free surface, a downhill particle flow is observed where 
the flow dynamics depends on the rotation speed of 
the cylinder 0,^. For low rotation speeds, the sur- 
face flow consists of individual avalanches called discrete 
avalanche regime. With increasing rotation speed the 
separation time of avalanches will decrease until no in- 
dividual avalanches are detectable. A nearly constant 
particle flow is found along the free surface and this 
regime is consequently termed continuous flow regime. 
Close to this transition, the surface can be well approx- 
imated by a straight plane, which can be used to deter- 
mine the surface angle. For even higher rotation speeds, 
cascading, cataracting and centrifuging particle motion is 
also observed |^]. Experiments performed in the discrete 
avalanche regime using mono-disperse particles show that 
the mixing time depends strongly on the volume filling 
fraction of the cylinder and no mixing is seen for an ex- 
actly half-filled cylinder A theoretical description 
could be given which is based on the mixing between 
wedges [^-||. On the other hand, when a mixture of 
particles which differ in size or density are placed in 
a rotating cylinder, the denser or smaller particles will 
concentrate in a central region close to the free surface 



after only a few rotations, which is termed radial segre- 
gation. This was studied experimentally and numerically 
for varying size ratios j^-p^ and density ratios [pl]-p^. 
The amount and direction of segregation depends on the 
rotation rate However, radial segregation is always 
observed in the continuous flow regime, regardless of the 
filling fraction of the cylinder and we are investigating 
numerically the interplay of the two competing phenom- 
ena of mixing and segregation, a problem first studied by 
Rose Q and restated by Behringer . We are partic- 
ularly interested in the dependence of the total amount 
of segregation on the size ratio of glass beads and how 
the segregation behaves for size ratios close to one. 

The paper is organized in the following way: After a 
brief motivation of the physical system in mind in the 
beginning of the next section, we will explain our numer- 
ical model and the physical meaning of the parameters in 
remainder of section ||. In order to model the dynamics 
of glass beads correctly with our model, we have to in- 
clude rolling friction. The details of our implementation 
are given in the appendix. We define a order parame- 
ter which allows to quantify the amount and the speed 
of radial segregation. The dynamics of this order pa- 
rameter are discussed in detail in section [II for different 



particle size ratios and as function of the volume filling 
fraction of the cylinder. In addition, we start with an ini- 
tial configuration similar to the one used in the original 
experiment on the mixing of mono-disperse particles |^ 
to demonstrate how the radial segregation competes with 
the mixing process. 



II. MOTIVATION AND NUMERICAL MODEL 

In many radial and axial segregation experiments, one 
of the components are glass beads. They are commer- 
cially available in large quantities and can easily be sieved 
to nearly uniform size distribution. When sufficiently 
large, the cohesion forces are negligible. In the contin- 
uous fiow regime the dynamic angle of repose for glass 
spheres is independent of the rotation speed of the cylin- 
der as long as the surface remains flat ||l^,|l3 . Since they 
are nearly perfectly round, particle rotations are an im- 
portant degree of freedom and cannot be neglected in 
a theoretical or numerical description. We use three- 
dimensional discrete element methods, also known as 
granular dynamics, this gives us the advantage to vary 
particle size ratios freely, whereas in experiments only a 
limited number of particle diameters and densities are 
available. 
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We want to investigate radial segregation, therefore we 
use only a cylinder with periodic boundary conditions 
along the rotational axis. In order to determine the min- 
imal length of this cylinder, to get no artifacts due to the 
periodic boundary conditions, we have investigated the 
range of the boundary effects in a previous study [0 . 



A. Forces during Collisions 

Each particle i is approximated by a sphere with radius 
Ri- Only contact forces during collisions are considered 
and the particles are allowed to rotate; we also include 
rolling resistance to our model to correctly describe the 
dynamics of glass beads (see appendix). The forces act- 
ing on particle i during a collision with particle j are 



F^j = -Y {Ri + Rj - rijfi) - ^nVijfi 
in the normal direction (ri) and 



Ff. 



-min(7,.% ■ s(t),^|i^"|) 



(1) 



(2) 



in the tangential direction (s) of shearing. In Eq. (|i|) 7„ 
represent the dynamic damping coefficient and Eq. (0) 
7s represent the dynamic friction force in the tangential 
direction, rij represents the vector joining both centers 
of mass, Vij represents the relative motion of the two 
particles, and Y is related to the Young Modulus of the 
investigated material. Dynamic friction in this model is 
defined to be proportional to the relative velocity of the 
particles in the tangential direction. 

During particle- wall contacts, the wall is treated as 
a particle with infinite mass and radius. In the normal 
direction, Eq. (^ is applied, whereas in the tangential 
direction, the static friction force 



min(fc, / V,, ■ sit)dt, ^i\Fl'\) 



(3) 



is used. This is motivated by the observation that when 
particles flow along the free surface, they dissipate most 
of their energy in collisions and can come to rest in voids 
left by other particles. This is not possible at the cylin- 
der walls. In order to avoid additional artificial particles 
at the walls we use a static friction law to avoid slip- 
ping and allowing for a static surface angle when the 
rotation is stopped. Both tangential forces are limited 
by the Coulomb criterion, see Eqs. (||) and (||), which 
states that the magnitude of the tangential force cannot 
exceed the magnitude of the normal force multiplied by 
the friction coefficient /i. Even though the experimen- 
tally measured friction coefficient for fresh glass beads is 
/i = 0.092 1 18 1 this can only be viewed as a lower bound 
in our case due to the wear of material caused by the 
uncountable bead collisions in the course of the experi- 
ment. For particle-particle collisions we use fj, = 0.19, 
and for particle-wall collisions, fiui = 0.6. The coefficient 
of restitution for wall collisions is set to 0.97 and to 0.831 



for particle-particle collisions, which are the measured 
values for glass beads |^8| and the density was set to 
p = 2.5^^. In order to save computer time, we set Y to 
6 • 10^ Pa, m which is about one order of magnitude softer 
than glass, but we checked that this has no effect on the 
investigated properties of the material . This gives a 
contact time during collisions of 1.1 ■ 10~^i, which is still 
quite small. The total number of particles we used were 
up to 17000. 



B. Order parameter 

To compare the quality and speed of the segregation 
of different runs with different parameters it is desirable 
to have only one order parameter to investigate all of 
these [Q. For this we divide the cylinder with Diame- 
ter D into n concentric hollow cylinders with thickness 
A — and measure the number density of the smaller 
particles, pi, in each hollow cylinder i ~ 1,... ,n (see 
Fig. |l|). We choose as such an order parameter q the sum 
over all positive deviations of pi above the mean number 
density po — (Number of small particles)/(/(i)/2)^7r) in 
a cylinder with a filling fraction /, normalized with re- 
spect to the ideally segregated case, where all inner cylin- 
ders i = 1, . . . , c are composed only out of small parti- 
cles shown in grey in Fig. |l] and in the outer cylinders 
i — c + 1, . . . , n only out of large particles. To get rea- 
sonable results, we take A to be of the order of the di- 
ameter of the larger particles, which gives in our case 
n = 1 1 . For smaller values of n the resolution gets worse 
and for higher n fluctuations become more pronounced, 
since there would be "empty" shells with no particle cen- 
ters in it. For varying n around 11, g has an over-all error 
of 0.03 for n e [9,20]. 
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1. Cross section through a more than half-filled cylin- 



To illustrate this procedure, we show in Fig. (where 
we used eleven concentric shells) the packing fraction of 
small (O) and large (-I-) particles in each cylinder for 
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a volume filling fraction of 66.4%. The volume filling 
fraction is defined as the the ratio of the volume occupied 
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FIG. 2. Volume fraction of small and large particles in each 
cylinder for a volume filling fraction of 66.4%. 

by the granular material to the cylinder volume, reading 



/ 



Voce 



, where the occupied volume depends on the 



packing fraction 77 via Vc 



-E 



Hence the packing fraction 77 is determined by 



all particles 3 



all particles 3 



i.e. filling the drum with a certain filling fraction and 
summing up the volumes of each particle. We found that 
this packing fraction rj does not change with the amount 
of achieved segregation. For Fig. ^ the particle size ratio 
of small and large particles $ := ^ was $ — 0-''5mm _ q g 

^ ^ R 1.5mm 

and the cylinder was rotated for 2^ revolutions which 
gives a nearly complete radial segregation. The small 
particles, denoted by O, are mostly found in the middle 
of the cylinder, whereas the large particles, denoted by 
+, show a higher concentration in the outer cylinders. 
Also shown is the total volume occupied by all spheres, 
denoted by □ , which gives an average value in each cylin- 
ders of 77 — 0.64 which is close to the value for a random 
particle packing. 



III. SIZE SEGREGATION 

For our simulation we are using a cylinder with diam- 
eter D — 7cm and periodic boundary conditions along 
the rotation axis. The cylinder length was 2.5cm which 
was sufficiently long to ensure that boundary effects are 



negligible p7| ] . The cylinder was filled with a binary mix- 
ture of large beads having a radius of i? 1.5mm and 
small beads r g {0.75mm, 1.0mm, 1.25mm}, where the 
small particles can have a concentration of 50% or 33% 
by volume. The aspect ratio of the drum diameter D 
and the average particle diameter 2r is D/ (2r) = 28, this 
is of the order of laboratory experiments, where we have 
D/{2r) = 25 up to 40 for the example in Ref. 0. The 
density of our beads is always p = 2.5-^^. For our value 
of f2 = 15rpm, we are in the continuous flow regime with 
a flat free surface and the Froude number is 

Fr = (2.3... 9.1) X 10-2 

so we can neglect inertia effects. 
1 
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FIG. 3. Typical time series of the order parameter q for 
two different filling fractions of the cylinder with 1.0mm and 
1.5mm beads. 



A. Time evolution of the order parameter 

Usually the initial state is a random mixture of small 
and large particles which gives a value of g « 0. For our 
value of f2 = 15rpm, we are in the continuous fiow regime 
with a fiat free surface and the order parameter will show 
a global trend of increasing in time and saturates on the 
long run when the cylinder rotation is started. A typical 
time evolution of q using a 50% volume fraction of 1.0mm 
smaller particles is shown in Fig. ^. The general trend 
can be well approximated by an exponential saturating 
function of the form 



a 



(4) 



with a characteristic segregation time tc and a final 
amount of segregation goo ■ The best fit to the data points 
for a 50% filling of each of l.Omni and 1.5mm beads 
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was obtained for the parameters tc = (6.1 ± 0.3)s and 
Qoc = 0.644 ± 0.040 which was added in Fig. |. 

When working in the discrete avalanche regime, it was 
found in experiments |^ and mathematical models [^|J|] 
that the least amount of geometrical mixing is given 
for a half-filled drum which could be explained by the 
avalanche mixing of wedges. Pcratt and Yorke |^ applied 
their model also to the continuous avalanche regime by 
taking the limit of an infinitely thin flowing layer and 
no change in the angle of repose during the flow. This 
should lead to "infinitely" many avalanches during one 
revolution and with a finite avalanche duration all these 
avalanches are superposed and could lead to continuous 
avalanches (steady flow). Motivated by these findings, 
we expected for our setup to get the least geometrical 
mixing for filling fractions around 50% as well and the 
strong fluctuations which are most pronounced for a fill- 
ing fraction of 50% are a clear sign for slow geometrical 
mixing. 



gives a = TT and l=2i? . From this we calculate (v) to 
be (v) — 12.5^ in each case. The depth of the fluidized 
layer is determined as in Ref. (2^j2l[] by looking at the 
velocity profile along a line through the center of the 
drum and perpendicular to the free surface. The depth 
of the fluidized layer is the distance from the free surface 
to the point where the velocity profile reaches its zero 
value. 

In Fig. ^, the thin line corresponds to a half-filled cylin- 
der and due to a slight asymmetric start configuration, 
which persists due to the bad geometric mixing, large 
fiuctuations which decrease in time are visible. For other 
filling ratios of the cylinder the asymmetry in the be- 
ginning gets erased by the geometrical mixing and the 
fiuctuations are suppressed or are decaying in a rapid 
way. This is illustrated by the thick line in Fig. [s] which 
is for a volume filling fraction of 80% and as expected the 
fiuctuations have a much smaller amplitude and do not 
show such a pronounced periodicity. 



1^1 




0.5 1 1.5 2 2.5 3 3.5 4 



T[s] 

FIG. 4. Fourier transform of q{t) of Fig. |^ . 

Fourier-transforming these fiuctuations gives main 
peaks for the 50% and the 80% filled drum at T = 
(2.56 ± 0.09)s and T = (3.18 ± 0.15)s, albeit the peak 
for the 80% filled case is much less pronounced than in 
the 50% case as is shown in Fig. ^; higher harmonics are 
also visible. These times are exactly the time it takes for 
a particle to make one revolution, i.e. to appear at the 
same spot again, given as 



where is the angular velocity of the drum, a the arc 
where the particle is in the solid block, / the length of 
the fiuidized layer and (u) the average velocity of the 
particles in the fluidized layer. For a half-filled drum this 



B. Dependence on the filling fraction 

We will now turn to the pre-factor of the exponential 
fit in Eq. (^, q^a, which quantifies the final amount of 
segregation and study its dependence on the volume fill- 
ing fraction and on the particle size ratio. This is shown 
in Fig. H for three different particle sizes, corresponding 
to a size ratio of <& = 0.5,0.67 and 0.83 and a concen- 
tration of small particles of 50%. For smaller particles, 
the final amount of segregation is higher for nearly all 
filling ratios of the cylinder. This can be explained by 
the higher mobility of the small particles which traverse 
through the network of voids of the large particles. 
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FIG. 5. Final amount of segregation for a concentration of 
50% of small particles and three different size ratios. 
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Since the geometrical mixing will cause also mixing of 
the segregated core of small particles with the large parti- 
cles, best segregation should be achieved for zero mixing, 
i.e. a filling fraction of 50% w hen mono-disperse particles 
are used, see also Sec. Ill A. In our case where we have 



a fluidized layer with finite width, the best segregation 
occurs not for a half-filled cylinder, instead it occurs for 
a cylinder where the solid block under the fluidized layer 
is 50%. The fluidized layer has a width of about three to 
four particles and hence the total filling fraction for the 
least geometrical mixing would be for a filling fraction 
of 60%, which is in agreement with our simulations (see 
Fig.|). 

Concerning different concentrations of small particles, 
our parameter was chosen in such a way that the fi- 
nal amount of segregation should be independent of the 
volume fraction of small particles in the cylinder. We 
checked this numerically for a concentration of 50% and 
33% and found a perfect agreement within the error bars. 
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FIG. 6. Characteristic number of revolutions for segrega- 
tion for a concentration of 50% of small particles. 

The speed of the segregation is characterized by tc, see 
Eq. (^, stating that for t = tc the system has reached a 
segregation of 63% of the final value of Qoo ■ In Fig. ^, we 
plot the characteristic number of revolutions, defined as 
ric = as a function of the volume filling fraction of 
the cylinder for three different sizes of the small particles 
corresponding to particle size ratios of $ = 0.50, 0.67 and 
0.83 again. The segregation times are sufficiently smaller 
for the smaller particles which is also found in experi- 
ments p2| . They also show a general trend of increasing 
with increasing filling fraction. However, for a more than 
half-filled drum, where the exact value depends on the 
particle size due to the different width of the fluidized 
layers, the segregation becomes faster again. This can be 
explained as follows: Due to the nature of the geometri- 



cal mixing, an unmixed core will persists in the middle of 
the cylinder for a filling fraction greater than 50%. The 
exact number depends on the width of the fluidized layer 
which depends e.g. on the rotation speed of the cylinder. 
Close to a filling fraction of one, only a small ring close 
to the cylinder wall can participate in the segregation 
process. Consequently, the final amount of segregation is 
small which agrees with Fig. |[ but this value is reached 
fast, therefore the segregation time, tc, is small. The 
numerical data indicates that the filling fraction which 
corresponds to the maximal value of tc decreases with 
decreasing particle size. 

It is well known, that the segregation process is faster 
and more pronounced if the particle size ratio becomes 
smaller ^2 . The results from a two dimensional rotating 
drum model indicate that segregation is observed for an 
arbitrary small size ratio , whereas the data from ver- 
tical shaking experiments suggest a cut-off ratio around 
$ = 0.5 1^^. In order to address this question, we show 
in Fig. 1^ the final amount of segregation, qoo , as a func- 
tion of the particle size ratio, <&, for a concentration of 
small particles and a volume filling fraction of 50%. Even 
though obtaining accurate data for values of $ close to 
one is rather difficult due to the long segregation time, 
the data shown in Fig. |^ support the hypothesis from |2|] 
that segregation will be present for any finite size differ- 
ence. This was determined by using a linear fit of the 
form 



goo($) = c(l-$o) 



(6) 



which gives c = 1.6 ± 0.1 and $o = 1-00 ± 0.02 when aU 
seven data points are used for the fit; shown as a dashed 
line in Fig. 
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FIG. 7. Final amount of segregation as function of the par- 
ticle size ratio ^ ~ (the vertical dashed line denotes "1>t). 

Due to our definition, the maximal achievable value for 
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the final amount of segregation is goo = 1- Therefore for 
small size ratios, the behavior must deviate from the lin- 
ear dependence which is already visible for the value for 
$ = 0.5 which was obtained by interpolating between two 
filling fractions. Obtaining data points for even lower val- 
ues of $ is nearly infeasible by todays computers, due to 
the large demand on computer time caused by the large 
particle numbers. For values of <f> < (wide dashed 

line in Fig. ^ we expect a completely different behav- 
ior since the small particles are then sufficiently small 
to propagate through the voids of a three-dimensional 
hexagonal packing. Please note that since we have a ran- 
dom packing, see Fig. the threshold value should even 
be higher and <f>T := "^^/^ J^^t serves as a lower bound. 



C. Half— filled pre— set cylinder 

After having discussed the segregation dynamics as 
function of filling fraction and particle size ratio in the 
preceding paragraphs, we will now turn to illustrate the 
interplay betwee n mixi ng and segregation. As already 
mentioned in Sec. Ill A, no geometrical mixing was found 
for an exactly half-filled drum in the discrete avalanche 
regime [^-|[. The continuous flow limit taken in Ref. ^] 
suggests that this is also true in the continuous flow 
regime but the numerical data sets given in Fig. || and ^ 
do not support this hypothesis which can be attributed 
to the fact that the fluidized layer has a non-zero width. 

In order to illustrate this point in more depth, we 
are starting with an initial configuration where the left 
half of the cylinder is purely composed out of large 
beads (R= 1.5mm) and the right half out of small beads 
(r= 1.0mm) giving a total number of 4420 particles (see 
Fig. H, top left picture) . After turning the drum counter- 
clockwise for 1.6 seconds at = 15rpm, which would sim- 
ple interchange the regions occupied by large and small 
particles if no mixing would be present, the interface is 
still well defined and nearly a straight line (top right pic- 
ture). After turning for 2.8 seconds, the interface be- 
tween the large and small beads is still quite sharp, albeit 
it is not a straight line anymore. After the start of the ro- 
tation, it takes 0.23 seconds for the continuous flow to set 



in and from Sec. Ill A, we recall that it takes roughly 2.65 
seconds for a particle to undergo a full revolution. The 
tongue of small particles into the large ones at the center 
of the cylinder is the starting point of forming a core of 
small particles. After rotating for 5.3 seconds, which cor- 
responds to two full particle revolutions, the formation of 
a core of small (white) particles is even more pronounced. 
After two full particle revolutions, t = 7.8s, the shape of 
the interface between large and small particles close to 
the cylinder wall becomes even more diffusive. It rather 
resembles a diffusion process along the azimuthal direc- 
tion which can be described in a similar fashion as the 
front propagation along the axial direction in rotating 
cylinders [p5[. The segregation mechanism in the pre- 



set cylinder starts immediately, i.e. no mixing of the two 
components is necessary in order to obtain a radially seg- 
regated core of small particles. The final picture in Fig. ^ 
(bottom right) corresponds to 28 particle revolutions and 
shows a nearly symmetric, well segregated cluster of small 
particles. Also note that hardly any large particles are 
found in the segregated core of small particles whereas 
smaller particles are still found close to the wall of the 
drum. We expect the latter effect to disappear when the 
drum is rotated for long enough times. 

In order to determine the degree of mixing in the hori- 
zontal cylinder along the initially sharp vertical front we 
use a procedure proposed by Metcalfe et al. ^ . We cal- 
culate the center of mass for each particle size, project 
it onto the free surface which is initially horizontal and 
calculate the distance of each of the two centers of mass. 




7.8 sec 56.2 sec 

FIG. 8. Different snapshots of the cylinder with a starting 
condition, where initially all the small (large) particles are on 
the right (left) side of the cylinder (particle radii 1.0mm and 
1.5mm). 

The time evolution of this distance, ^c, which was made 
dimensionless by dividing by the distance of the start 
configuration is shown in Fig. ^ It corresponds to the 
configuration shown in Fig. Ig, i.e. a half-filled cylinder 
containing an equal volume fraction of 1mm and 1.5mm 
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particles. In Ref. |3|, this procedure was used to show 
the mixing of mono-disperse particles in a rotating drum 
at a filling fraction of / = 39%. Since geometrical mixing 
is observed for this filling fraction, the centroid positions 
decayed in time and could be well approximated by 

CeW=cos(^)e^*/*= (7) 

where T stands for the period, see Eq. (^. On the 
contrary, no geometrical mixing is observed for mono- 
disperse particles in a half-filled drum which would cor- 
respond to a characteristic time tc — oo in Eq. (0) and 
lead to non-decaying oscillations. For a binary particle 
mixture, the segregation process will lead to a decay of 
the distance of the two centroids in time for any filling 
fraction and we have chosen to present numerical results 
for counter-clockwise rotation and / = 50% in Fig. || to 
illustrate this. 
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FIG. 9. Normalized centroid position projected onto the 
free surface of a half- filled cylinder with 1.0mm and 1.5mm 
particles. 

The numerical data can be well fitted by an exponen- 
tially decaying oscillation according to Eq. . This gives 
T = 2.48s which is in excellent agreement with Fig. ^, and 
values of tc — 20.7s for clockwise rotation and tc = 20.1s 
for counter-clockwise rotation which are the same within 
the error. The exponentially decaying part was added 
to Fig. ^ as dashed line. In the beginning, t < 2Gs a 
slight asymmetry is visible towards negative values which 
was also observed in other numerical simulations using a 
two dimensional geometrical model However, if this 
persists on the long run and will led to a final non-zero 
value for the centroid position could not be determined 
by this procedure since the deviations from the exponen- 
tial fit were usually higher than the calculated offset of 
(0.004 ± 0.001)cm. Since the numerical data seems to 



indicate that a pure exponential decay is too slow in the 
beginning and two fast for longer times we also tried a 
stretched exponential decay of the form e~'*/*'='''^ . This 
gives values of tc = (14.0 ± 1.8)s and (3 = 0.88 ± 0.05. 
However, we can not rule out that the deviation of the 
exponential law may be due to the numerical noise and 
a more general statistical theory is needed to resolve this 
question. 

The procedure described above to follow the centroid 
dynamics is not capable to determine the depth of the 
centroid position below the free surface due to the pro- 
jection onto the surface. However, the last picture of 
Fig. H shows that the number of layers of large particles 
below and above the segregated cluster of small parti- 
cles is not the same. This leads to different distances 
of the center of mass below the free surface for small 
and large particles. To determine its dynamic, we plot 
in Fig. |l^ the Euclidean distance of the two centroids. 
When properly shifted by T/2, the curves for clockwise 
and counter-clockwise rotations show a similar behavior 
where fluctuations are a little less pronounced for the lat- 
ter case. The minima correspond to configurations simi- 
lar to the one shown in the bottom right picture of Fig. ^ 
when the two centers of mass lie on a line which goes 
through the origin of the drum. Even though the mini- 
mal distance during the time evolution is close to zero, a 
non-zero value of 0.22cm is estimated for the stationary 
state without oscillations. 
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FIG. 10. Distance of the centers of a half-filled cylin- 
der with 1.0mm and 1.5mm particles ( - - - clockwise, — 
counter-clockwise rotation) . 

IV. CONCLUSIONS 

The main conclusions of this paper are as follows. In 
the continuous flow regime size segregation takes place for 
arbitrary small differences in particle size. To quantify 
this result we introduced an appropriate order parameter. 
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which allowed us to compare directly all different drum 
scenarios. From this we showed, that the radial segrega- 
tion process is faster and more pronounced for particles 
with a large size difference and that the final amount of 
segregation shows a linear dependence when approaching 
a particle size ratio of one; thus no threshold value for ra- 
dial segregation exists, which was an unresolved question 
for long time. 

We also studied in detail the interplay between mix- 
ing and segregation in rotating cylinders for different 
volume filling fractions of the cylinder and found that 
the highest achievable segregation can be obtained for a 
slightly more than half-filled cylinder and therefore also 
least mixing. The difference to a mono-disperse system 
in the discrete avalanche regime could be attributed to 
the fluidized layer leading to a partial destruction of the 
underlying already segregated core where the destruction 
increases with layer width. 

When starting with an initial configuration that con- 
tains well-separated regions of small and large particles 
no mixing of the components is necessary in order to ob- 
tain a radially segregated core. This inter-penetration 
process resembles a diffusion process and segregation 
starts immediately without undergoing a previous mixing 
of the two particle components. 
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APPENDIX: ROLLING FRICTION 

One of the most distinct properties of a glass bead is 
its ability to roll, therefore rolling had to be included 
into our model. The drawback was, that a glass bead 
would have had a Coulomb friction of zero with our fric- 
tional laws; i.e. a glass sphere would start to roll even 
on an infinitesimal inclined plane or a particle would roll 
on forever on every flat plane, which in reality clearly 
does not occur. For an ideal elastic particle, the defor- 
mation of the particle would be symmetric to the point 
of contact and therefore the resulting counter force of 
the plane would be exactly opposite to the gravitational 
force for all times. In reality the deformation is not 
elastic, i.e. the deformation lags behind as is indicated 
in Fig. ^ for a particle on an ideal hard plane. The 
counter force of the plane F acting on the particle gets 
mediated by the deformation of the particle. The point 
where this force acts on is shifted slightly by ro to the 
back, the normal component of F compensates the grav- 
itational force exactly (otherwise the particle would 



bounce); leaving the tangential component of F which 
acts as rolling resistance which must be compensated by 
a dragging force for a particle with constant velocity on a 
flat plane. Also simulation shows that the angle of repose 
is much too small in comparison with experiment with- 
out rolling resistance. To overcome this weakness we add 
rolling resistance to our model, see Fig. Ill], by using 



R 



F, 



N 



(8) 



Here the rolling resistance ro is a constant material pa- 
rameter and results from the slight viscoelasticity of the 
materials, rg is of the order of 10~^ — 10~^mm for most 
the materials. For particle-particle interactions, we take 
the same law for rolling resistance as for particle- wall in- 
teractions. 



R 



\ Fn. 




^Fn J 









FIG. 11. Viscoelastic rolling sphere on a hard surface. 

The rolling resistance acts as a net torque constructed 
out of a force couple Fr with 



F. 



iF^I^(nxs). 



(9) 



We also have to consider that the rolling resistance can 
only decrease angular momentum, but never revert it. 
And so we have to limit F^ by the quantity that would 
reduce the angular momentum to zero within the next 
time step, namely 



'-mR ■ {{h X s) ■ Lo)/{M). 



For the torque we therefore get 

T — R ■ min(_FV, -Fr, 



J- 



(10) 



(11) 



As can already be seen from the definition of the rolling 
resistance, Eq. (H), the rolling friction on small particles 
will be higher than on large particles which was also ob- 
served by studying one particle on a bumpy line and 
in experiments with glass marbles. To illustrate this fact, 
we show in Fig. the dependency of the angle of repose 
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on the rolling friction parameter tq. The higher friction 
of smaller particles results in a steeper slope of the angle 
of repose in Fig. An important result of this is that 
we can adjust the rolling friction in such a way that small 
and large glass beads have the same angle of repose as 
is seen in experiments for glass beads Also one 

clearly sees that for small tq the slope of the angle of re- 
pose 1^ is proportional to tq , which result from the law 
of rolling resistance Eq. (||). 




0.001 0.002 0.003 0.004 0.005 0.006 

ro [cm] 

FIG. 12. Angle of repose for different friction coefficients 
and particle diameter of 1.0mm and 1.5mm for three different 
values of /I (/I = 0.1, 0.2 and 0.3 from bottom to top). 
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